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ABSTRACT 

The mean velocity and Reynolds stress tensor throughout a whirling annular seal are presented. The data was collected with 
a three dimensional laser Doppler velocimeler using phase averaging. Two axial flow conditions (Re = 12,000 and 24,000) were studied 
at one shaft speed (Ta = 6,600). The eccentricity and whirl ratios were 50% and 100%, respectively. There is a region of high axial 
momentum at the inlet on the pressure side of the clearance that migrates around the seal to the suction side at the exit. The normalized 
axial momentum in this region is higher in the low Reynolds number case due to an axial recirculation zone that occurs on the suction 
side of the rotor at the inlet. The recirculation zone docs not occur in the high Reynolds number case. At both Reynolds numbers there 
is a recirculation zone on the rotor surface in the pressure side of the inlet. This recirculation zone extends from 20 to 200* past the 
rotor zenith in the tangential direction, and is one third of a clearance wide radially. The high Reynolds number recirculation zone is 
1.5 mean clearances long, while the low Reynolds number zone extends 2 mean clearances downstream. When compared to previous 
studies, it is apparent that the tangential momentum is no greater for a seal with whirl than for one without if other parameters are 
constant. Areas of high tangential momentum occur in the clearance where the axial momentum is low. Average exit plane tangential 
velocities in the high Reynolds number case are 1.5 times greater than those in the other flow case. These results are in general 
agreement with predictions made by other investigators. 


INTRODUCTION 

Annular seals control leakage from high to low pressure areas in pumps, compressors and other turbomachines operating at 
high speeds. A secondary but equally important purpose is to provide rotordynamic stability (Tam et al., 1988). The seal separates high 
pressure from low pressure regions, but some leakage is desirable, providing damping for the vibrating system, cooling on the shaft, and 
increasing stability. Seals with a small axial Reynolds number (generally below 2000) are generally unstable while seals with large axial 
Reynolds numbers produce large damping and stiffness coefficients and therefore contribute to stability (Allaire et at., 1978). 

The mechanical performance of an annular seal is implicitly dependent on the fluid flow through the clearance, since the forces 
in the seal are generated hydrodynamically. Destabilizing forces tend to move the rotor of an annular seal from its centered position, 
and whirl usually follows. Whirl occurs when the center of the rotor precesses around the center of the stator. The whirl ratio is defined 
as w=u)Ju) l . The rotor excitation can become so great that the rotor will contact the stator. 

There is usually much effort devoted to characterizing and predicting the damping, stiffness and inertia variables in a seal to 
decide whether it will be stable prior to production. While the macroscopic effects have been studied extensively, the flow field within 
the seal has not been investigated because of the difficulty of the measurement. The clearances must be very small to simulate a real 
machine, and it is very difficult to measure velocities in the wake of the rotor. 

Previous Research Annular seals have been investigated since the mid 1970’s, when they were identified as contributors to turbomachinc 
instabilities. In 1977 data became available for a wide range of operating conditions because of research directed toward the fuel pumps 
in the Space Shuttle’s main engine (Tam, 1988). The following is a review of some relevant research regarding fluid flow through eccentric 
annular seals. The investigations mainly involved rotordynamic modeling, but one effort was directed wholly to measuring the flow field 
in an annular seal. 

Lessen (1987) performed an analytical study of the flow in a dynamically eccentric whirling annular seal at an arbitrary Taylor 
number within the Taylor vortex regime. Axial flow through the seal was modeled as plug flow which was justified by assuming high 
momentum transport from Taylor vortices. Boundary layers near the rotor and stator surface were dominated by Goertler disturbances 
(instabilities akin to Tollmein-Schlichting waves). Goertler disturbances were induced by the curvature of the boundary surface, and in 
this ceometry are of far less importance than the Taylor vortices (Schlichting, 1979). The width of the Goertler disturbance-dominated 

regime was taken to be one unit of the non-dimensional Blasius variable, at which point the local velocity u!U tvz was 0.33. Lessen 

suegests that marginally unstable Goertler disturbances will be superseded by Tollmein-Schlichting waves if the sum of squared axial and 
tangential Reynolds numbers exceeds the critical Reynolds number for their formation: Re,'4-Re # "> Re CTlt . This indicates that the nature 
of the turbulence in the seal resembles flat plate turbulence al high tangential and axial Reynolds numbers. It was further found that 
the reaction forces in the radial direction were inertial in nature and could induce instabilities, and that the tangential reaction force was 
purely a function of viscous forces. Finally the ratio of potential core circulation to shaft circulation was the same in the eccentric rotor 
case as in the concentric. Lessen’s important finding was that tangential flow circulation did not change when the rotor was offset 
eccentrically (whirl was not considered). 


1 This work was performed under NAG3-181 with the NASA Lewis Research Center and the Texas A&M University Turbomachinery 
Laboratory Research Consortium. 



Chen and Jackson (1984) studied the effect of axial eccentricity (displacement in the r -0 plane) and misalignment (angular 
displacement in the Z-r plane) on rotordynamic coefficients in annular seals. They suggested that the forces generated in high pressure 
seals are greatest in straight annular seals, and developed relationships between reaction forces and leakage rate through the seals that 
included the effects of eccentricity, misalignment and rotation. A concentric tapered annular seal was used as a model for the eccentric 
seal, and it was found that the effect of eccentricity or misalignment on the leakage rate was not as great when the flow regime was 
turbulent. The seal leakage increased with the degree of eccentricity, and decreased with the degree of misalignment. 

Allaire et al. (1978) obtained an analytical solution of the semi-empirical bulk flow equations developed by Hirs (1978) for 
pressure in an eccentric non-rotating annular seal. The seal length-to-diameter ratio was 0.16 and the axial Reynolds number was about 
an order of magnitude larger than the tangential. He found that in seals with a small axial Reynolds number the pressure drop though 
the seal was mainly due to friction effects, but as the leakage rate increased the Bernoulli effect became stronger and high pressure 
gradients existed through the seal. At very high flow rates, large head losses occurred at the entrance to the seal and the pressure 
distribution throughout the seal tended to be constant. The load-carrying capacity of the seal was small for the low flow rate case 
increased as the now rate increased, peaked where the pressure gradients were high, and then decreased as the flow rale increased 

Hashimoto et al. (1988) investigated short, eccentric journal bearings to find the effect of fluid inertia on trajectory This was 
basically an investigation of whether laminar flow theory would suffice in turbulent conditions. The length-to-diameter ratio of the journal 
was 0.5. The averaged form of the momentum and continuity equations was used in the analysis and reaction forces were calculated for 
the fluctuating pressure term. Through numerical simulation the trajectory of the journals was calculated and it was demonstrated that 
fluid inertial effects help contribute to the stability of the system. 

Kanemori and Iwatsubo (1989) considered a long annular seal with a length-to-diameter ratio of 3.0. The eccentricity of the 
shaft was adjustable and misalignment was not considered. They concluded that for small axial flow rates the results were similar to 
predictions from Fritzs theory. This did not consider the axial flow and simply suggested that the reaction force in the seal was 
dependent on the whirl velocity. It was determined that the tangential reaction force was destabilizing for whirl ratios between zero and 
1/2, and that the seal was unstable at small Re x , but stable at higher Reynolds numbers. The amplitude of the tangential reaction force 
was minimized when the whirl ratio was 1/2. 

Simon and Frene (1989) studied convergent and divergent annular seals with a cryogenic working fluid. They considered the 
compressibility and variable viscosity of the fluid, inlet swirl and pressure drop effects. The initial assumption was that inertial effects 
of the fluid in the seal could be neglected. A numerical example was run with axial Reynolds numbers about 130,000 and 160,000 This 
example indicated the utility of the variation of the fluid properties in the analysis. Their results were in good agreement with 
experimental data. ^ 

Tam et al. (1988) conducted a numerical study of eccentric annular seals without misalignment, using a model of the fluid forces 
based on the average fluid circumferential velocity ratio. The model was based on the assumption that the dynamic forces were rotating 
at the precession speed. The numerical grid was three dimensional and rotated with respect to a stationary observer at the precession 
rate of the rotor, i.e. the rotor was a stationary reference frame. The grid consisted of 12x6x16 cells in the tangential, radial, and axial 
directions, and was considered invariant in time since the orbit was circular and periodic. The Prandtl mixing length turbulence model 
was selected because of its simplicity and because the mixing length was obviously the clearance of the seal. The eddy viscosity concept 
p,-p/V/>, was used to simulate the effect of the turbulence levels on the mean velocities. The pressure drops ranged from 0.3 to 18.5 
bars, the rotational shaft speeds from 1442 to 5085 rpm, the eccentricity ratios from 0.24 to 0.8, and the whirl ratios from 0.3 to 1.0. 
Bromotriflouromeihane and oil were modeled as the working fluids. 

The results of the study indicated that there were significant changes in the local values of the seal dynamic forces and that 
arge tangential separation zones existed through the seal. The existence of the recirculation zones was found to depend on the shaft 
and precession speeds, with increasing likelihood at lower shaft speeds and higher precession speeds (whirl ratio approaching 1.0). 
Recirculation zones were located along the stator wall, beginning just prior to the minimum clearance location, and rotated at the 
precession speed. Negative preswirl (opposing the direction of seal rotation) intensified the recirculation zones while preswirl in the 
direction of rotation weakened the secondary flow. 

The effects of eccentricity and rotation speed on leakage flow rate were also studied. Leakage was most dependent on 
rotational speed when the precession speed was low. Increasing precession speed increased the leakage rale through the seal because 
of a pumping effect. 

Fluid injection from four stator positions into the clearance was also studied. It was found that injection in the direction of 
the shaft rotation increased leakage while injection against the rotation decreased leakage. Injection essentially mimicked the preswirl 
effects mentioned earlier, but downstream of the inlet. 

Johnson (1989) and Morrison et al. (1991) used a three dimensional laser Doppler anemometer to study an annular seal 
positioned first concentrically and then eccentrically. The eccentricity was created by moving the stator closer to the rotor at an angular 
position of interest (the rotor was not whirling). Johnson measured axial-radial profiles at 0=0*, 90°, 180°, and 270°. The leakage flow 
was^4.86 1/s (77 gpm, Re = 24,000) and the rotational speed varied from 0 to 3600 rpm(Ta = 0 to 6,600). The preswiri was 0", -45*, and 

For the rotor concentric with the stator, tangential velocities were higher (~0.3 W >fc ) for positive preswirl than for no preswirl 
(-0.18W*) or negative preswirl (-0.25W* at the inlet and 0.10VV b at the exit). For the eccentrically mounted stator, tangential velocities 
vary to satisfy continuity considerations being greater in the small clearance than in the large. In the large clearance the tangential 
velocity did not increase along the length of the seal but in the small clearance there was azimuthal acceleration. The effect of preswirl 
was similar to that for the concentric geometry. 

In summary, previous studies showed that eccentrically aligned rotors affect annular seals in several ways. The core circulation 
in an eccentric seal is the same as in concentric seals (Lessen, 1987). Whirling rotors increase the leakage flow rate through seals 
presumably due to a pumping effect (Chen and Jackson, 1984 and Tam, 1988). High leakage rales through the seal decrease the axial 
pressure gradient across the rotor because of high head losses in the entrance region (Allaire, 1978). High turbulence levels enhance the 
stability of whirling journals, suggesting that once whirling begins it may be damped out (Hashimoto et al„ 1988). Seals are unstable al 
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low axial Reynolds numbers (below 2000) (Allaire et al. 1978 and Kanemori andlwalsubo, 1989). Reaction forces are destabilizing when 
the whirl ratio is less than 0.5 (Kanemori and Iwatsubo, 1989). Tangential recirculation zones discovered in whirling seals during certain 
operating conditions were seen to diminish if there was preswirl or fluid injection in the direction of rotor rotation, and were enhanced 
if preswirl or injection opposed the rotor circulation direction (Tam, 1988). Morrison, et al. (1991) and Johnson (1989) showed that 
preswirl effects tangential velocity development in concentric and statically eccentric seals. 

Measuring the velocity profile in the whirling seal is the only means to test the validity of assumptions used in the analysis of 
other researchers. For example, it is highly unlikely that the axial velocity distribution in the seal can be modeled as a plug flow as 
hypothesized by Lessen. Theories that were proposed in these studies and tested include the proposition that 1) the tangential velocity 
distribution does not increase in the presence of whirl, 2) there is a tangential velocity recirculation zone along the stator, and 3) that 
seal leakage increases with eccentricity. 


PROGRAM OBJECTIVES 

To address the lack of experimental data, the mean velocity profile and turbulence levels in an annular seal (37.3 mm long, 
164.1 mm diameter, and a nominal clearance of 1.27 mm) were measured for evaluation of flow phenomena. Figure 1 illustrates the test 
section in which water is the working fluid. The stator contains an optical window through which the beams of a three-dimensional laser 
Doppler velocimeter are directed. An eccentricity ratio of 0.5 is imposed on the rotor by mounting it on an eccentric bushing, which is 
attached to the shaft of an electric motor (the eccentricity ratio is defined as the ratio of the off-centeredness of the rotor to the clearance 
width, e = e/c). A detailed description of the facility and error analysis are contained in Johnson (1989) and Morrison, et al.(1991). In 
real machines eccentricities can be attributed to misalignment during assembly, side loads such as rotor weight, impeller loads (Nelson 
and Nguyen, 1988), or destabilizing seal forces. Instabilities caused by pressure fluctuations in the flow can also induce eccentric whirling. 

By varying the flow rate through the seals, Reynolds numbers and Taylor numbers typical of full scale applications were 
obtained. The seal angular speed was 3600 rpm for all the tests. The leakage flow rates were 2.43 1/s (38.5 gpm) and 4.86 l/s (77 gpm), 
with an eccentricity ratio of 50%. Since the rotor was connected directly to the shaft of the motor, the whirl ratio was 100%. The 
Reynolds numbers investigated were 12,000 (U mv ,=3.7 m/s) for the first flow case and 24,000 (U 1V ,=7.4 m/s) for the second, while the 
Taylor number for both cases was 6,600 (W lh = 28.7 m/s). 


RESULTS 

The rotor is spinning clockwise and whirling in the same direction as it is rotating. Angles are measured from the minimum 
clearance point in a clockwise direction, and the minimum clearance point is the rotor’s zenith. The suction side of the clearance is on 
the left of the rotor, and the pressure side is on the right. 

Re= 12.000 

Mean Velocities Significant features of the mean velocity profiles are that the axial velocity peaks on the pressure side of the clearance 
at the seal entrance and rotates around the seal circumferentially until it is on the suction side at the exit. The radial velocity distribution 
generally points outward to the stator at the seal inlet, and dissipates downstream. An axial vena contracta exists over the first two 
clearances on the pressure side of the rotor. 

Figures 2 and 3 present contours for the three mean velocity components and the axial velocity variance. Figure 2 represents 
the flow at five axial locations along the entire length of the seal while Figure 3 presents data from -0.02 to 0.14 Z/L of the flow. 
Upstream of the inlet (Figure 3) the flow is very sensitive to the motion of the rotor, as indicated by the large radial velocities away from 
the rotor where it accelerates toward the stator. On the pressure side of the seal there are large tangential velocities and very small radial 
velocities. The large radial velocity zone is explained partially by the flow in the plenum upstream of the seal moving with the rotor from 
top to bottom. The rotor is pulling the flow along with it for angles larger than 180° and pushing flow away from it for angles less than 
180°. Once the flow enters the seal clearance, the high radial velocities at 90° diminish and the flow follows the movements of the rotor 
more closely. At the inlet the dfd9 gradients are most apparent. No tangential recirculation zones were seen along the stator in this flow 
case; therefore, if one exists it must be between the stator and the first radial grid line. 

Upstream of the inlet to the seal, at Z/L=-0.02 (Figure 3), the maximum axial velocity (1.3U #W where U tv .=3.7 m/s) occurs +95° 
from the rotor zenith (0 = 0° position). There is a negative zone of axial velocity on the other side of the seal which reaches -0.4U 1W (1.5 
m/s). It is evident that at this location the rotor is affecting even the flow upstream of the seal. The radial flow component is moving 
very quickly away from the rotor (0.9U lve -3.3 in/s) from about 0=90° to 180®, but is depressed everywhere else. This is representative 
of the motion of the rotor, which is pulling fluid toward it in the lower half of the stator, and pushing it away on the other side. The 
tangential flow component is uniformly distributed around the clearance at 0.12W lk (W lk is the rotor surface velocity, 28.7 m/s), which 
is on the order of U. The inlet preswirl is caused by the rotation of the rotor which influences the upstream velocities through the large 
volume of water trapped in the hollow section of the rotor. The tangential velocity is accelerated into the minimum clearance area 
beginning at 260°. To satisfy continuity, the axial or radial velocities must be reduced. Consequently, the negative axial zone begins at 
260° also, but the radial velocities do not seem to react. 

At the entrance, Z/L=0 (Figures 2 and 3), the maximum axial velocity increases to 2.4U lve (8.9 m/s) and the peak location 
moves to 70°. The size of the back flow area diminishes, but it does not move relative to its previous position. It appears that the peak 
axial velocity is located in a region of decelerating tangential flow. For equal axial pressure drop across the seal at all azimuthal angles, 
the greatest axial velocity should be in the largest clearance area if there is no rotation or whirl. But, the tangential component 
accelerates toward the minimum clearance area, so the axial momentum decreases to satisfy continuity. The smallest tangential velocities 



occur in a band from 60° to 230°, where the largest axial velocities are located. The greatest radial velocities continue to occur between 
90* and 180°, and the maximum radial velocity is still 0.9U iVt (3.3 m/s) outward. 

The axial velocity component accelerates down the seal, and the peak begins to rotate further away from the rotor zenith at 
Z/L=0.04 (Figure 3). The maximum axial velocity reaches 2.6U 1W (9.6 m/s) and rotates to 90° as the tangential momentum transports 
the large axial momentum zone around the seal. At the same time, an axial recirculation zone (due to the vena contracta effect at the 
seal entrance) extends around the rotor, from 80 to o20° on the rotor surface. The shear layer between the positive and negative axial 
velocities develops intensely at 130*, and the large axial velocities progress outward from the rotor toward the stator side of the clearance. 
The tangential velocities begin to increase also, especially on the suction side of the seal, where there is little axial momentum. The radial 
velocities clearly define the vena contracta region on the pressure side of the clearance, a strong counter flow exists at the rotor surface 
(«350° to 180°) with a positive radial velocity further into the clearance. The large positive radial velocity zone has moved around the 
seal slightly to 300 s . There are large dV/dr gradients in the shear layer region around the axial recirculation zone and through the vena 
contracta shear layers. 13 

Farther down the seal, at Z/L=0.07 (Figure 3), the U component begins to drop in magnitude, with a peak velocity of 2.1U lvt 
(7.8 m/s) at 0=80°. The recirculation at 180° on the bottom of the rotor decreases in magnitude, and the recirculation on the opposite 
side of the rotor from the maximum velocity region subsides as well. The large U zone begins to spread out around the seal, and is 
located out closer to the stator than the rotor. Large «V stresses along the rotor surface divert the flow out to the stator. The radial 
velocities are directed into the rotor in much of the zone occupied by the high axial velocities, but are positive in the U recirculation areas. 
The W distribution is more uniform, at 0.20W lk (5.7 m/s) throughout much of the seal. The tangential velocity contours are closer 
together in the small clearance area of the seal, but they are still evenly distributed throughout the clearance. It appears that the dW/dr 
gradients scale with the clearance width. 

Between Z/L=0.07 and 0.21 the peak axial velocity decreases and spreads over a much greater area of the clearance as the 
flow continues downstream. The deceleration is due to shear along the stator wall and conservation of momentum, since the high axial 
now area is spread out over more of the seal. This is the region where the vena contracta diverges out to reattach to the rotor surface, 
and the axial recirculation disappears. The steep dVfd r gradients toward the rotor have also diminished at this point. The positive and 
negative radial velocity zones begin to equalize on either side of the seal and approach zero. 

At the midplane of the seal (Z/L=0.50, Figure 2) the U profile becomes almost uniformly distributed, with a peak of 1.1U #VB 
(4 m/s) at 0 = 135". The 0.8U JW , (3 m/s) contour extends throughout most of the circumference of the seal. The radial velocities are 
virtually zero throughout the entire seal, and the tangential velocities remain very uniformly distributed at an average of 0.20W >k (5.7 m/s). 

Down the seal at Z/L=0.77 a radical shift in the axial velocities occurs. The peak velocity increases to 1.5U #ve (5.5 m/s) and 
rotates around to 0 = 295*. The axial recirculation does not reappear, and the peak velocity is fairly close to the rotor, indicating that there 
are low turbulence levels there. The positive and negative radial velocity zones have also rotated to opposite sides of the seal, with the 
negative areas following the peak axial velocity zone. The magnitude of the radial velocities is below 0.05U 1¥t (0.18 m/s). The tangential 
profile shifts with the larger tangential velocities located in the pressure side of the seal, which is consistent with the downstream trend 
that large tangential velocities belie small axial velocities, and vice versa. 

The axial velocity profile continues to change all the way to the exit plane (Z/L= 1.00) with the peak axial velocity progressing 
toward the rotor zenith. The radial velocities are similar at Z/L = 0.77 and 1,00. At the exit, tangential velocities reach a maximum 
across the width of the clearance with high tangential velocities extending from the pressure side of the seal at 0 = 60* to the suction side 
at 0 = 200° at values near 0.3W >h (8.6 m/s or 2.3U tw J. 

Turbulence levels Upstream of the inlet to the seal the u ' m ' Reynolds stress is evenly distributed throughout the clearance averaging 
0-1 U,v. (Figure 3, 1.4 m7s fc ). This stress level is what remains of the grid turbulence generated at the entrance to the plenum. At the 
inlet, the normal stresses begin building up, thereby inhibiting the flow along the rotor. The u ' u ’ stress increases sharply along the rotor 
from 45° to 175°, and the du T u ' /dr gradients are very large halfway into the clearance. This normal stress, which is analogous to a 
pressure (Rodi, 1984), retards the axial flow along the rotor. Similarly, gradients of the normal stress are considered as dynamic pressure 
gradients (Tennekes and Lumley, 1990). The high "pressure" area blocks axial flow along the rotor, diverts it up into the vena contracta, 
and induces recirculation behind it. Downstream of the inlet at Z/L=0.04 (Figure 3), the u ' u ' stresses reach a maximum of 0.9U 2 (12 
m7s 2 )at 70°. The peak is located a quarter clearance from the rotor, begins developing at 31CT, builds to a maximum at 70*. then 
diminishes to a minimum at 300* . The du' u ’ /dr gradients between the u ’ u ' peak and the U peak are very steep in the shear layers 
along the recirculation zone. Three quarters of the clearance from the rotor and at 80* the mu’ stresses recede to 0.2U, W 2 , and the U 
velocity reaches a maximum. This location is where the vena contracta reaches its narrowest width. The du ' u ' id$ gradients are greater 
on the upstream side of the peak and then trail off at a gentler slope downstream. At Z/L=0.07 the u ' u ' stress is convected tangentially 
around the sea! to 80* from 70° just upstream. Th e struc ture is still intact but the gradients are smaller than before. The peak is still 
dropping, reaching 0.7U‘ (9.5 nr/s 2 ) at 80* The u ’ u ' stresses dissipate rapidly between Z/L=0.07 and Z/L=0.21. The peak level 
convected to 100 , so the convection of the u li stress increases also. The u " >_ m ' gradients in both the axial and tangential directions 
decrease in magnitude. At Z/L=0.21 the u u normal stress is distributed across the entire clearance on the pressure side of the rotor. 
This reduces the mean velocity peak, whereas before the high m ’ m ’ levels on the rotor surface had pushed the high U velocities toward 
the stator. At the midplane of the seal (Z/L=0.50) the u ’ u ' stresses stop changing any further. The peak u 1 m ' stress is 0.15U tvt 2 (2 
nr/s“), and it is located in a bubble on the bottom of the rotor where the recirculation was located nearer to the inlet. Another bubble 
develops in the minimum clearance area. For the rest of the seal these structures change little, suggesting that production and dissipation 
are in equilibrium from this point on. 

Re =24.000 


Mean Velocities The mean velocity contour plots reveal an interesting change in the upstream flow patterns with the increased How rate 
through the seal. Upstream and at the inlet to the seal (Figures 4 and 5) the flow vectors on the pressure side are pointing inward, as 


opposed to outward as they did in the slower flow rate. UpstreamofthesealatZ/L=-0.02(Figure5),themeanaxialveIocityprofiIe 

reaches a peak of greater than 0.9U 1V€ (6.6 m/s, where U„ e = 7.4 m/s) between 60° and 135°. There is no backflow out of the seal as in 
the lower Reynolds number case because the upstream pressure is higher, and the pressure developed in the seal cannot overcome it. 
The radial velocities are generally pointed inward toward the seal as the flow has just been diverted inward by the curve in the plenum 
wall. The greatest inward mean radial velocities are located in the same region as the high U’s. There is much less significant preswirl 
at the inlet for this case. The inlet W does not exceed 0.08W lh (0.33U 1V „ 2.4 m/s) over most of the clearance. 

At the seal inlet (Figure 4) the axial flow component begins accelerating into the seal. The peak axial velocity is 1.5U ave (11 
m/s) across a band from 30° to 60* and stays above 0.5U lve (3.7 m/s) across the entire circumference of the seal. The V profiles continue 
the strong inward trend. The maximum inward velocity reaches -0.8U >vt (-6 m/s) from 45* to 190*. It appears that the inward momentum 
imparted to the flow by the curved stator wall is still affecting the fluid, but at the same time the inlet shape of the seal is directing flow 
to the cavity under the seal. The tangential velocities begin to increase, with the 0.08W th (2.4 m/s) contour extending almost to the stator 
from 280° to 60°, and moving outward through the suction portion of the clearance. 

Just downstream of the inlet at Z/L=0.04 (Figure 5) the axial velocities are much greater, with a maximum of 2.2U iW (16.3 m/s) 
occurring at 95°. This peak is very close to the stator, and there are very high dU/dr gradients on either side of the peak. A shear layer 
exists from 90* to 180* and stretches away from the rotor on the maximum clearance side. A mild recirculation zone exists between the 
rotor and the shear layer with a maximum backflow velocity of 0.1U #vt (-0.7 m/s). This recirculation zone occupies less than 15% of the 
clearance and is not visible on the figures. Reattachment of the axial velocity occurs at 180* where the axial velocities are significantly 
smaller than on the other side of the clearance. The axial recirculation produces a small vena contracta similar to that for the low 
Reynolds number case but much smaller in magnitude. The V profile changes significantly at this point, as the direction of the velocities 
is redirected back out toward the stator. A maximum of 0.1SU Jve (1.3 m/s) at 45* is surrounded by mild dV/d r gradients on either side. 
The radial velocities are pointed outward near the rotor and inward near the stator on the wide clearance side of the rotor, indicating 
the presence of a vena contracta around 135*. Because the axial momentum in this area is so small, the tangential velocities are able to 
develop quickly, especially in the recirculation zone where W = 0.24VV ih (6.9 m/s, 0.9U Bve ) in the middle of the clearance from 60* to 280*. 

From Z/L=0.04 to Z/L=0.07 (Figure 5) the axial velocities decrease, only reaching a maximum of 1.9U #ve (14 m/s) at 0=120*. 
The high velocity area is still along the stator, and the steep OV/ot gradients near the rotor surface are relaxing. The recirculation zone 
has disappeared along the rotor, but the velocities are still very' low there. The disappearance of the shear layer implies that the axial 
streamlines arc reattached to the rotor. The radial velocities are directed back out toward the stator after turning at the reattachment 
point. The tangential velocity development diminished with the reattachmcnt of the axial velocity streamlines. 

Moving to Z/L = 0.14 (Figure 5), the U profile is more evenly distributed throughout more of the seal, as the peak axial 
momentum region reaches wider clearances, which is decelerating the flow. The peak axial velocity is 1.5U 1V . (11 m/s) and it is located 
at 90*. The rest of the profile remains attached to the rotor and stator with very mild d\J/dr gradients. The axial velocity on the suction 
side of the rotor is steadily increasing, with the minimum continuous contour at 0.7U <ve , where just upstream it is 0.6U iVt . The radial 
velocities continue to decrease in magnitude and now arc fairly insignificant. The decrease in axial momentum is acting favorably for 
the tangential velocities, which are now increasing over the entire seal. The 0.08W lh (2.4 m/s) level is now extending much further into 
the maximum clearance region, and it is distributed more uniformly. 

From Z/L=0.14 to Z/L=0.50 (Figure 4) the trends discussed above continue; the axial velocity continues to spread more 
uniformly across the entire seal, the radial velocities remain small, and the tangential velocities increase with the decrease in axial 
momentum. At the midpiane the axial profile is almost uniform, with regions in the maximum clearance region having velocities which 
are a little higher. As the location of the maximum axial velocity rotates around the seal, the negative radial velocity zone also moves. 
By the midplane it is relocated to the minimum clearance region. ITie tangential velocity develops to 0.08W tll all the way out to the stator, 
and it is increasing in the pressure section of the clearance, reaching 0.l6W th (4.8 m/s) by the midplane. 

Farther down the seal at Z/L=0.77 the location of the maximum axial velocity is convected around the seal by the increasing 
tangential velocity component. The peak axial velocity is at 240“ and has increased to I-4U ave (10.4 m/s) as the tangential component is 
forcing the fluid into a smaller region. The radial velocity is similar to that upstream, with positive values pointing out from the rotor 
and negative in from the stator. The tangential profile develops significantly, with the 0.24W lh (6.9 m/s) zone extending almost entirely 
across the pressure zone around 60*. This significant increase in the tangential momentum is the mechanism driving the axial velocity 
peak around to the other side of the clearance. 

The axial velocities on the suction side of the clearance continue to increase as the flow progresses downstream while the 
tangential velocities increase on the pressure side. This continues to the exit plane where the peak axial velocity reaches 1.7U #Vf . (12.6 
m/s) at 265*. At the exit plane nearly the whole pressure side of the clearance is dominated by tangential velocities of at least 0.24W jh 
(0.93U) and on the suction side they are 0.16W jh (4.6 m/s). The radial profile remains virtually unchanged through the rest of the seal, 
not exceeding 0.08U lvt (0.6 m/s). The radial velocities are pointed inward along the stator. 

Turbulence levels Upstream of the inlet, in the plenum at Z/L = -0.02 (Figure 5), the normal Reynolds stress term is extremely subdued, 
which is reflected in the mild gradients in the mean velocity terms. Hie u r u ' term barely reaches more than 0.05U 1V / (2.7 m 2 /s 2 ). At 
the inlet (Figure 4) the from edge of the rotor causes changes. 4Tie u u term builds (exceeding 0.33U #W< \ 18 nr/s 2 )along the rotor 
surface on the bottom side where the recirculation zone is located a bit farther downstream. (Note: the recirculation is evident in the 
numerical data but does not appear in the contour plots.) This stress decelerates the axial velocity at the surface of the rotor. By 
Z/L=0.07 (Figure 5), the u ' u ' value decreases with the large gradients also showing a marked reduction. The peak of the u ’ u y ~ stress 
is transported to 110* by the tangential momentum, and the peak has diminished. Dissipation dominates the flow at this point, where 
upstream turbulence production is the major mechanism, and works to distribute the turbulence term fairly evenly by midplane. The 
profile remains reasonably uniform on the suction side, but by ihe exit plane. Z/L=1.00, there are again high values of u u Although 
the gradients are steep, it seems likely that the phenomenon is real because of the circumference over which it is distributed. The 
turbulence producing source at this plane is the step at the end of the seal. 


SUMMARY 


The now field inside an annular seal whirling at a whirl ratio of 1 and an eccentricity ratio of 0.50 has been measured using 
a 3-D laser Doppler anemometer system. Reynolds numbers of 12,000 and 24,000 were investigated while the Taylor number was held 
constant at 6,600. The flow structure, as evidenced by the velocity and turbulence distributions in the whirling annular seal proved to 
be far different than in a statically eccentric seal otherwise operating at the same conditions. Prominent features included a peak axial 
velocity that began in the pressure section of the clearance at the inlet and rotated around to the suction side of the clearance at the exit, 
a vena contracta on the pressure side of the rotor, and no significant increase in tangential velocity with whirling rotor motion when 
compared to statically eccentric seals. 

The axial flow profiles for both the high and low Reynolds number cases have similarities but also display substantial 
differences. In both the profiles, the peak axial velocity rotates from the pressure to the suction side of the clearance. In the low 
Reynolds number case, this started at 60” with a magnitude of 2.4U lvl over the inlet. By midplane the profile spreads uniformly across 
most of the clearance. This rotates around the clearance, and at Z/L=0.77 the peak reaches the final circumferential position of 300”. 
Downstream of this position, the peak continues to grow, increasing to 1.9U TO at the exit of the seal. 

The mean velocity field reveals a highly three dimensional flow with large radial velocities near the inlet of the seal. The 
normalized axial momentum near the inlet on the pressure side of the rotor is higher in the low Reynolds number case due to an axial 
recirculation zone that occurs on the suction side of the rotor at the inlet. Another recirculation zone exists on the rotor surface on the 
pressure side of the inlet at both Reynolds numbers. This recirculation zone extends from 20“ to 200“ in the tangential direction, and 
is one third of a clearance wide radially (not apparent in the contour plots). The high Reynolds number recirculation zone is 1.5 mean 
clearances long, while the low Reynolds number zone extends 2 mean clearances downstream. When compared to previous non-whirling 
flow studies (Johnson (1989)) at similar flow parameters, the tangential momentum does not significantly increase with the onset of whirl 
Areas of high tangential momentum occur in regions of the clearance where the axial momentum is low. Average exit plane tangential 
velocities in the low Reynolds number case are 1.5 times greater than those in the high Reynolds number case. 

NOMENCLATURE 


A Leakage area, nr 

c Mean clearance between the stator and rotor, m 

d Rotor diameter, m 

e Eccentricity, m 

L Length of the seal, m 

\ m Mixing length, m 

Q Leakage flow rate, m 3 /s 

r Radial distance from stator centerline, in 

Re Reynolds number = 2Uc/v 

Re *> Reynolds number in axial (x) or tangential ( 0 ) direction 

Re cr.i Critical Reynolds number for transition to turbulence 
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Axial mean velocity, in/s 
Bulk mean velocity=Q/A, m/s 
Radial mean velocity, m/s 
Tangential mean velocity, m/s 
Tune averaged Reynolds stress, nr/s“ 

Rotor surface speed, in/s 

Distance downstream of the seal inlet, m 

Eccentricity ratio, e/c 

Azimuthal angle measured in the direction of rotation frot 
Fluid absolute viscosity, (N s)/nr 
Turbulent eddy viscosity, (N s)/m : 

Strain rate 

Kinematic viscosity, nr/s 
Fluid density, kg/m 3 
Whirl ratio, u/^at T 
Rotor precession speed, rpm 
Rotor speed, rpm 


the minimum clearance 
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